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We discuss a supersymmetric unified model based on a product gauge group SU(5) xSU(5) xSU(5), 
where the gauge symmetry breaking is achieved without the adjoint or higher-dimensional Higgs 
field, and the doublet-triplet splitting in the Higgs masses is realized by the use of the discrete 
symmetry. In this article we present an explicit model for realistic fermion masses with the discrete 
symmetries Z7 x Z2- It is shown that all the observed masses and mixing angles for quarks and 
leptons, including neutrinos, are well described by the breaking of the symmetries imposed in the 
model. Especially, the maximal and large mixing angles in the atmospheric and solar neutrino 
oscillations are obtained as the most preferred values, and the typical value of the neutrino mixing 
element Ue3 is 0.1-0.3. We also point out the non-trivial relations among the /^-parameter for the 
Higgs mass, the charged fermion hierarchies, and the neutrino masses. These relations suggest that 
the scale of fi is of order of the weak scale. 



PACS numbers: 12.10.-g, 12.15.Ff, 14.60.Pq 

I. INTRODUCTION 

The unification of strong, weak and electromagnetic 
interactions is a very attractive idea ,1,]. The gauge 
coupling unification supports strongly this grand unified 
theories (GUTs) with a supersymmetry at low energy. 
The simplest model is based on an SU(5) group 2]. In 
spite of the theoretical appeal, there have been several 
puzzles in construction of supersymmetric GUT models: 
The first puzzle is the Higgs sector. In general, the very 
complicated and innovative Higgs sector should be con- 
sidered. The GUT gauge breaking usually requires the 
Higgs field(s) in the adjoint or higher-dimensional repre- 
sentation(s). In particular, the most annoying difhculty 
is the mass splitting between Higgs weak-doublets and 
color-triplets. Various mechanisms for this problem have 
been proposed [1 0,1111 0- 

The second one concerns the matter sector. Quarks 
and leptons are unified into irreducible representation(s) 
of the GUT group. This matter unification ensures the 
electric charge quantization, but also predicts a promis- 
ing mass relation, rrib = m-r, in the simplest SU(5) model. 
However, it predicts also the unwanted mass relations, 
uis = rUf^ and = mg, which are badly broken in na- 
ture. Moreover, recent neutrino experiments 8. Q. llollTH 
reveal that the neutrinos own non-zero masses which are 
much smaller than the electron mass, and remarkably 
that the mixing angles in the atmospheric and solar neu- 
trino oscillations are both large. The rest angle is only 
limited by experiment ,1^. We then have to explain 
small and large mixing angles in the quark-lepton sec- 
tor at the same time. The observed fermion properties 
may also require the GUT models to be complicated. 

One way to evade these issues is to give up the funda- 
mental assumption - a simple group -, i.e., to go into the 
product group like GxG 01 GxGx G, where G denotes 
the GUT gauge groups [IlEl El El 111 111 El 111 111. 
As pointed out in Ref. [ll| the class of these models. 



called the product-group unification, still preserves the 
nice features of the GUTs. Namely, the standard gauge 
group is contained in the diagonal subgroup, and hence 
the coupling unification is realized at tree-level. In addi- 
tion, quarks and leptons are unified in each GUT group 
G, and the electric charge quantization is also explained. 

There are three major advantages in the product-group 
unification: First, the breaking of gauge symmetry is ac- 
complished without introducing the Higgs fields in the 
adjoint or higher-dimensional representations |l3j |. This 
seems very promising for the connection to the string 
theory. In fact, it is known that there is no state in the 
adjoint or higher representation in the string spectrum 
with the afhne level one . Second, the doublet-triplet 
problem may receive a natural solution by the use of 
the discrete symmetry El E3 El El ES-^ Finally, we 
may gain important insights to the flavor physics. It is 
beyond the scope of the usual GUTs (as well as the stan- 
dard model) to understand the origin of the hierarchical 
structure in the Yukawa couplings, which are essentially 
free parameters. On the other hand, the mass hierar- 
chy between different flavors may be explained by the 
breaking of the product group when they belong to dif- 
ferent GUT group factors [lj|- This illustrates the idea 
of Refs. |3j l25l | in the unified theories. Moreover, the 
discrete symmetry for the doublet-triplet splitting can 
be utilized as the horizontal flavor symmetry, which is 
also crucial for the suppression of the dimension-five op- 
erators to avoid a rapid proton decay 0, 0, El |^ . 

In this article, we present an explicit model for 
realistic fermion masses based on the product group 
SU(5)xSU(5)xSU(5). Especially, we emphasize on the 
neutrino properties. Although the similar analyses have 



^ In the different class of models with a product group the unbro- 
ken R-symmetry is used for the doublet-triplet splitting l23l . 
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been done in Refs. O [Tl|, the model in Ref. Q 
leads to the small angles for neutrino mixing, and in 
Ref. 19] the neutrino physics has not been discussed. 
Our model is constructed by imposing the discrete sym- 
metries Z7 X Z2 , where the Z-j symmetry is introduced to 
solve the doublet-triplet splitting. The observed masses 
and flavor mixing of quarks and leptons, including neu- 
trinos, are well described by the breaking of the product 
gauge group together with these discrete symmetries. It 
is shown that the Z-j breaking is crucial for generating 
(i) an effective /i-term for the Higgs doublets, (m) the 
charged fermion mass hierarchies, and (iii) the Majorana 
masses for right-handed neutrinos.^ We point out that 
the /i-parameter at the weak scale is suggested from not 
only the observed masses of charged fermions, but also 
the neutrino masses indicated by the oscillation experi- 
ments. 

This article is organized as follows: in the next sec- 
tion, we first present our model based on a product group 
SU(5)xSU(5)xSU(5). In particular, the gauge symme- 
try breaking and the doublet-triplet splitting are dis- 
cussed. In Sec. mil we construct the fermion mass ma- 
trices from the invariance under the gauge symmetry as 
well as the discrete symmetries. The masses and mixing 
angles of quarks and leptons, including neutrinos, are 
discussed in Sec. IIVI We compare the prediction of our 
model with the experimental data. In Sec.[V]the impact 
on the leptogenesis of our model is considered. Finally, 
Sec. IVII contains our conclusion. 



II. SU(5)ixSU(5)2xSU(5)3 MODEL 
A. Gauge Symmetry Breaking 

In this article we will investigate a supersymmetric 
model based on a product group: 

G = SU(5)i X SU(5)2 X SU(5)3. (1) 

One of the advantages in this class of models based on 
product groups is that the gauge symmetry breaking can 
be achieved without introducing Higgs fields in the ad- 
joint or higher-dimensional representations. In fact, as 
pointed out in Ref. ^3|> the bifundamental fields are 
sufficient for this task. We introduce here three vector- 
like pairs of the bifundamental fields, (Ti,Ti), (Di,Z?i) 
and (V2,V^2)- Their gauge charges are presented in Ta- 
ble^ Although Ti(Ti) and Di(Di) have the same gauge 
charges, but they carry different charges under discrete 
symmetries. The breaking of the group G into the stan- 
dard model group SU(3)xSU(2)xU(l) is achieved by the 



TABLE I: Field contents of the model. 
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vacuum expectation values (VEVs) of these fields: 

(Ti) - (Ti) ^ diag( Ti , Ti , Ti , , ) , (2a) 
(i?i> - (:Di> -diag(0, 0, 0, (2b) 
(I/2) = (Fa) =cliag(T2, T2, T2, D2, (2c) 

where we take all the VEVs are real and positive."^ 
These bifundamental fields play also an important role 
for fermion masses and mixing. As we will see in Sec. IIVI 
the observational data suggest the scale of the VEVs is 
typically of the order of 0.1 (Af* is the cut-off scale 
which is taken as the reduced Planck scale in this anal- 
ysis). Below these scales, we obtain the minimal super- 
symmetric standard model. 

The above VEVs might be determined by a superpo- 
tential for the bifundamental fields. Such a superpoten- 
tial should avoid the appearance of the massless Nambu- 
Goldstone fields. The standard model group is embedded 
in a diagonal SU(5) group of G, and the gauge coupling 
unification is ensured at the leading order, although we 
have to include the threshold corrections around the G 
breaking scale. In this article, we will not discuss these 
issues, but we assume the VEVs as in Eqs. (0. 



B. Doublet- Triplet Splitting 

One of the serious difficulties in construction of GUT 
models is the doublet-triplet splitting in Higgs masses. 



^ The relation between the large Majorana masses for right-handed 
neutrinos and the breaking of the horizontal global symmetry has ^ We use sometimes the same letters for the superfields and for 
been discussed in Ref. |2bI. their VEVs. 
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There exist two weak-doublets of Higgs fields, Hu and 
Hj^, in the minimal supersymmetric standard model, 
where Hu and Hd are superfields whose VEVs give 
masses to up- and down-type quarks, respectively. We 
incorporate these Higgs doublets by introducing 

H = iGu,H^f, H^{G'd,Hdf, (3) 

together with color-triplet Higgs Gu and GJ^. Here H 
and H are 5 and 5*-plets of SU(5)3 and SU(5)i, respec- 
tively. It should be emphasized that the gauge symmetry 
G forbids an invariant mass term for H and H, however, 
the terms HTiH and HDiH are allowed in the superpo- 
tential. Although the former term is desirable to give a 
heavy mass to the color-triplet Higgs, the latter should 
be forbidden to leave the weak-doublet Higgs massless 
below the G breaking scale. How can we distinguish the 
weak-doublet Higgs from the color-triplet Higgs in a nat- 
ural way? 

The doublet-triplet splitting can be realized by using 
a discrete symmetry in the framework of the product- 
group unification [H Q El El H ISll. Here we shall 
follow the discussion of Ref . [21| . Such a symmetry can 
be constructed from a usual discrete symmetry Zj^ and 
a discrete "hypercharge" subgroup of U(l)i C SU(5)i or 
11(1)3 C 811(5)3. In the present model we choose that of 
U(l)i for example. Its element (combining with the Zj^ 
transformation) is parametrized as 

1- / -1 -1 -1 HiS. ii+l N / ,x 
gi = diag( a , a , a , , ) , (4) 

where a is an A^-th root of unity. For the later purpose 
we consider only an odd integer N. We assign the Z^ 
charges for the bifundamental fields as 

1 I ^ a 1 1 , Ui ^ a 2 Ox, 

and V2 and V2 carry zero charge. Then, it turns out 
that the VEVs shown in Eq. Q preserve a Zj^ symmetry 
which is the combination of Z^ and gi. 

This unbroken discrete symmetry Z^r can realize the 
doublet-triplet splitting. We assign the Z^ charges for 
H and i?: 

H^a''"H, H^a~''»^^H, (6) 

with an integer jih (= 0, • • • , A — 1), and the transfor- 
mations under the Z^ of each components are given by 

We can see that the Z^ symmetry, which is left after 
the G breaking, discriminates between the weak-doublet 
Hd and the color-triplet GJ^ as long as A 7^ 5. In fact, 
if it is the case, the term HDiH will be absent in the 
superpotential, while the following term is still allowed 

W ^HTiH , (8) 



which generates a heavy mass of Ti ~ 0.1 to the color- 
triplet Higgs. 

The exact Zn symmetry leads to no mass term for the 
weak-doublet Higgs in the superpotential. The so-called 
/i-term which is needed for the electroweak symmetry 
breaking is generated only by its breaking. For this pur- 
pose we introduce a gauge-singict field $ which carries 
a Zat charge -1-1 and includes the following term (hence- 
forth we take the cut-off scale Af* to be one) 

W ^ c^TlDiH^^ , (9) 

where is a constant and we set = 1 for simplicity. 
The effective /u-parameter then becomes'' 

JV + 5 

/i = Di$— . (10) 

For examples, when /i = 100 GeV and Di = 0.1, the 
VEV of $ is given by $ = 1.4 x IQ-* (A^ = 3), 2.7 x 10"^ 
{N = 7) and 6.3 x lO^^ (N = 9). 

C. Matter Fields 

Now we turn to discuss the matter fields in our model. 
In the simplest SU(5) model one family of quarks and lep- 
tons are grouped into a 10-plet and a 5*-plet. Similarly, 
we introduce the matter fields as 10-plets and 5*-plets 
of an SU(5) factor of G, which guarantees the electric 
charge quantization. Concretely speaking, 10^, 5* and 
li [i — 1,2,3) are introduced as the matter fields. The 
charges under the gauge group G can be also found in Ta- 
ble |I| First, each 10^ is allocated to a 10-plet of SU(5)i, 
respectively. Then, the gauge charges of H and H results 
in a unique choice for the distribution of 5* to cancel the 
gauge anomalies. Namely, one of 5*, say 5^, should be 
a 5*-plet of SU(5)2 and the rest, 82 and 83, should be 
5*-plets of SU(5)3. Finally, three gauge-singlet fields li 
(right-handed neutrinos) are introduced to generate neu- 
trino masses. 

This matter content leads to interesting features of 
the model: First, apart from the Majorana masses for 
right-handed neutrinos, the gauge symmetry G allows a 
Yukawa coupling only for top-quark. All other Yukawa 
couplings are induced effectively by the G breaking and 
are suppressed by some powers of VEVs of the bifunda- 
mental fields. This accounts for, at the first approxima- 
tion, a large top-quark mass and smaller masses for the 
other fermions. 

Second, our model offers naturally the so-called "lop- 
sided family structure" mill m [33. The distribution 
of lOi under G is different from that of 5* to cancel the 
gauge anomalies. This may be the reason why there is a 



^ Although the effects of the supersymmetry breaking may also 
give an additional contribution to /x, it can be neglected when 
the supersymmetry breaking field carries the Zjv charge zero. 
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difference in the mass hierarchies between up- and down- 
type quarks. More importantly, the same gauge charge 
between 82 and Sg may ensure a large v^-Vt mixing angle 
confirmed in the atmospheric neutrino oscillation.^ 



III. FERMION MASS MATRICES 

In this section we present fermion mass matrices in our 
model. As discussed above, the product gauge group G is 
considered as the flavor symmetry. The discrete symme- 
try for the doublet-triplet splitting is also important 
for fermion masses. Since the bifundamental fields, H 
and H carry non-trivial Zj^ charges, the matter fields 
should also transform non-trivially under the Zjv. Some 
Yukawa couplings are forbidden and induced effectively 
by its breaking together with some power of the VEV of 
the field This means that the symmetry is crucial 
not only for the doublet-triplet splitting but also for the 
hierarchies in fermion masses. 

In this article we present the model based on the Zj 
symmetry (iV = 7). This is because, first, the doublet- 
triplet splitting is realized only if iV ^ 5. Second, when 
N = 3, the VEV of $ is too small to explain the charged 
fermion masses [see the discussion below Eq. (I1U|| ]. Fi- 
nally, the larger TV, i.e., the larger VEV of $, is disfa- 
vored since it induces too small neutrino masses. In the 
following, therefore, we will construct the fermion mass 
matrices based on (i) the gauge symmetry G and (ii) the 
Zj symmetry for the doublet-triplet splitting. Further, 
we will impose (iii) an additional Z2 symmetry. The 
reason for its inclusion will be clear later. 



A. Mass Matrices for Charged Fermions 

Let us first consider the Yukawa terms for charged 
fermions, and we will discuss neutrino masses later. We 
define the Yukawa couplings by 

W = H uluj + Y,HjHdldj + Ye^jHe'^e, , (11) 

where Ui and are the left- and right-handed up-type 
quarks, etc. Without taking into account the discrete 
symmetries, the gauge charges of matter fields give the 



effective Yukawa couplings: 



Y, 



Y, = 



Y 



Di Ti D2 T2 1 

T1T2 Dip2_ D1T2 

Ti D1D2T2 Di 
Ti D1D2T2 Di 

D1P2 D1D2 D1D2 
Di Di D2^ Di 
Di Di D2^ Di 



(12a) 



(12b) 



(12c) 



where the bifundamental fields should be considered as 
their diagonal components which have non-zero VEVs 
(cf. {T2) = (T2) = T2 and {D2) = (D2) = D2). These 
results are obtained at the leading order of the power 
expansion in the bifundamental fields. As we will study 
in Sec. lIVI all VEVs of the bifundamental fields turn out 
to be of order 0.1 (in the unit A/* = 1), and hence the 
leading terms are sufficient for our discussion. 

It should be noted that we have omitted dimensionless 
coefficients in front of each element in the above Yukawa 
couplings. In principle, we can take them as any com- 
plex numbers. Here, however, following to the idea of the 
Froggatt- Nielsen mechanism Isi), we assume that these 
unknown coefficients are of order one, and thus the hi- 
erarchy in the Yukawa couplings are generated by the 
breaking of the product gauge group G (as well as the 
discrete symmetries, see below). Therefore, we should 
keep in mind that our analysis is of order of magnitude 
wise and receives the uncertainty from our ignorance to 
determine these coefficients. 

As mentioned above, the Z-j symmetry for the doublet- 
triplet splitting gives rise an additional hierarchical struc- 
ture in the Yukawa couplings H12|l . Due to the smallness 
of the suppression factor, <I> ~ 10"'^, we have to choose 
carefully the Z7 charges for matter fields as well as Higgs 
fields in order to obtain realistic fermion mass matrices. 
Notice that the Z7 charges for bifundamental fields have 
already determined in Sec. ^ 

Let us start with the charge assignment of the third 
family fermions. For a large top-quark mass we require 
the following term in the superpotential without any sup- 
pression factor of $ 



T4^ D IO3 IO3 , 
and also we need the term. 



W D HDi IO35 



3 ' 



(13) 



(14) 



^ The same structure of the matter and Higgs fields (but the dif- 
feren t st ructure of the bifundamental fields) can be found in 
Ref. I19I . However, the different mass matrices are induced and 
also the implication to the neutrino properties are not discussed 
there. 



to explain bottom and tau masses. Moreover, the neu- 
trino mass would come from the term, W D H H 5^5^, 
which, unfortunately, induces too small neutrino mass 
scale compared with the neutrino oscillation experiments. 
From this reason, we require the following term as a con- 
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sequence of the seesaw mechanism |3a.l3^:^ 



W D -HH5pl. 



(15) 



With these three requirements we may choose the Zn 
charges as 



H -.2, H -.A, IO3 : 6, : 2 . 



(16) 



The charges for other matter fields are found as fol- 
lows: We find from and nis that the charges for IO2 
and 62 should be same as IO3 and Sg, respectively. Of 
course, we have to avoid too small electron mass when 
the charges of lOi and 5^; are determined. The results 
are 



10i:l, 102:6, 5^:3, 5^:2. 



(17) 



These charges are also listed in Table |3 With this Z7 
charge assignment the Yukawa couplings take the forms: 



/ 



Y, = 



Y 



\ ^DiTi D2T2 1 

T1T2 D1D2 D1T2 
$Ti D1D2T2 Di 
$Ti D1D2T2 Di 

f $DiD2 <^>^DiD2 <i>^DiD2 

'i'^Di DiD2^ Di 
\ ^^Di DiD2^ Di 



(18a) 



(18b) 



(18c) 



Comparing with Eqs. H12() there appear suppression fac- 
tors in some Yukawa couplings. 

We find, however, the Yukawa couplings in Eqs. H18|) 
predict wrong mass ratios, mu/mt, mc/mt and md/mh, 
which conflict with the observation by orders of magni- 
tude. The present model calls for an additional discrete 
symmetry Z2 to correct these mass relations. We assign 
non-trivial charges only for the bifundamental and Higgs 
fields. Since all matter fields carry even parity of Z2, 
this Z2 symmetry should not be considered as a usual 
horizontal flavor symmetry. Notice that this Z2 symme- 
try does not affect the mass term for the color-triplet 
Higgs in Eq. ((SJ and the effective /i-term in Eq. 
From the Z2 symmetry the parity-odd components in 
the Yukawa couplings receive a suppression factor of 
(DiDi) = Tt{DiDi) = 2Df. In the following analysis 
we neglect this factor of two which is beyond our approx- 
imation. 



The same assignment of the Z7 charges might induce the opera- 
tor W D H H 5^5 ^, h owever, our model is not such a case as 
we will show in Eq. I47i . 



Finally, the effective Yukawa couplings in the model 
are given in the forms 



Y,. 



Y. = 



Y,. 



I _^Di{p^Di)_ T^D2T2 tA 

^ D,_T,T2HD,Di) D2{p^_pi) T2A , (19a) 
\ (^DiTi{D^D^) D2T2 I I 

fTiT2{Di^i) $62^1^2 ^'^i'lTsX 

^Ti{D{Di) D1D2T2 Di , (19b) 
\^Ti{DiDi) D1D2T2 ) 

( ^T3iD2 Di'D2{DiDi) <^'^ DiD2{D{Di)\ 

<^'^Di{DiDi) DiD2^ Dx 
\<S>^DiiDiDi) DiD2^ Di ) 

(19c) 



By putting the non-zero VEVs, the matrices for charged 
fermions are obtained as follows: 



Me' = Hd 



( $L>i3 Ti^D2T2 Ti^ 
^Di^TiT2^ Di^D2 T2^ 
\ ^Di^Ti D2T2 1 

Di^TiT2 ^^Di^D2 ^^Di^T2 
$Ti D1D2T2 Di 
$ri DXD2T2 Dx 

^DiD2 <i>^DiD2 ^^DiT2 
^^Di^ DiD2^ Di 

$21)^3 DiD2^ Di 



(20a) 



(20b) 



(20c) 



It is rather important to mention that the second and 
third rows in Md and M^^ have the same structure, since 
there is no symmetry which distinguishes between 62 and 
83. In this analysis, however, we assume that unknown 
coefficients in front of 62 and 63 are different. Note that 
some coefficients are related of each other due to SU(5) 
gauge groups. For example, the 3-3 components in Md 
and Mf, are identical. 



B. Neutrino Mass Matrices 

Now, we turn to discuss mass matrices for neutrinos. 
We define the Dirac and the Majorana mass matrices by 

W = vl {Mnh ^2 + \ {MNh ^1 ■ (21) 

Since right-handed neutrinos = li are totally gauge 
singlets, the gauge group G tells nothing about the struc- 
ture of their Majorana masses. On the other hand, they 
can carry non-trivial charges under the discrete Zj sym- 
metry for the doublet-triplet splitting, and the origin of 
the Majorana masses and their structure may be ex- 
plained by the Z^ breaking. Moreover, this fact leads 
to the relation between the /z-parameter for the Higgs 
mass and the neutrino masses in our model. As we will 
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show in the next section, this relation fits very well to 
the observations. 

We assign here the Z-j charges for the right-handed 
neutrinos in the following way: 



li : 2, I2 : 3, I3 : 3. 
Then, the Majorana mass matrix becomes: 



Mat = $ 



/ $2 $ $ 
$ 1 1 

U 1 1 



(22) 



(23) 



The Z-j breaking generates the Majorana masses Mi (i 
1,2,3) as 



Ml 



M2.3 



(24) 



This shows that right-handed neutrinos acquire super- 
heavy masses due to $ = C(10~'^), hence they offer a 
natural setup for the seesaw mechanism . 

The Z'j charges for right-handed neutrinos in Eq. I|22|) 
are also important to construct the Dirac mass matrix 
My). These charges, together with the gauge symmetry 
G and the discrete symmetry Z2, generates in the 
form 



( Di^D2 $ $ 
<^>^Di^D2 1 1 
y <^>^Di^D2 1 1 



(25) 



Here the suppression factor Di^ appears in the first row 
due to the Z2 symmetry. Similar to the mass matrices for 
down-type quarks and charged leptons in Eqs. if^ : the 
second and the third rows in Md take the same structure 
because 83 and Sg carry the same charges under all the 
symmetries. These facts lead to the large flavor mixing 
between the second and third families in the left-handed 
leptons and right-handed down-quarks. Furthermore, we 
should note that the 2-1 and 3-1 elements in Md are ex- 
tremely suppressed and essentially zeros. These missing 
elements give us the desired hierarchy in the neutrino 
mixing angles. We will give these particulars in the next 
section. 

To summarize this section, we have constructed the 
mass matrices for all quarks and leptons including neu- 
trinos based on the product group G, the Zi symmetry 
for the doublet-triplet splitting, and the additional Z2 
symmetry. The structure of the mass matrices arise from 
the breaking of these symmetries under our assumption 
that all the unknown coupling constants be of order one. 
In the following section, we will compare the obtained 
results with the observation. 



IV. FERMION MASSES AND MIXING 

We are now at the point to evaluate masses and mixing 
angles of quarks and leptons in our model. The fermion 



mass matrices are shown in Eqs. (|20|) . (|23|) and H25|l . 
Apart from unknown coefficients of order one, they are 
all determined by the following six parameters: 



Ti, Di, T2, D2, tan/3. 



(26) 



where tan/3 — Hu/Hd, thus we expect rather non-trivial 
predictions between masses and mixing angles. 



A. Charged Fermions 

We first discuss masses and fiavor mixing of charged 
fermions. From Eqs. (|20|l the third family masses are 
found to be 

mt ^ Hu, mi,~mr ^ Ha Di , (27) 

and we also find the following mass hierarchies: 



:::i^D^D2, ^-02X2, ^^02', (28a) 
nit m-b m-T 



mt nib 



Di Ti T2 , 



$1)2. (28b) 



Note that masses in the first family can be estimated in 
a rather solid way since the off-diagonal elements in the 
first column of M„ and in the first row of Md and MfP" 
are strongly suppressed. On the other hand, the elements 
of the CKM matrix are found as 

n nV ' \^^b\^D2T2. ^ . (29) 

Ul U2 ±2 J^l 

Here \Vus\ and \Vub\ come almost from the mixing in 
down- type quarks. 

It should be noted that the above relations between 
masses and mixing angles are valid at the unification 
scale. Of course, we should also keep in mind that these 
relations are predicted with uncertainty coming from un- 
known coefficients in the mass matrices. 

As shown in Eq. (|27|l , the bottom-tau unification mb = 
rrir in the simplest SU(5) model holds within a good ac- 
curacy even though the mass matrices Md and M^^ in 
Eqs. H2U|) are different. In addition, the unwanted SU(5) 
mass relations, = rus and rrie = md, are avoided in 
the model, although quarks and leptons are unified into 
irreducible SU(5) representations. This is because of the 
gauge charges of matter fields as well as the structure of 
our vacuum For instance, the Georgi-Jarlskog rela- 
tion (34I, rUfj, ~ 3771s, can be realized by taking D2 — 3T2. 
Moreover, we have a factorization of the mixing angles 
for quarks. 



WusWVci 



I Kb 



(30) 



since the quark mixing originates mostly in the structure 
of the mass matrix for down- type quarks. 

All the VEVs of bifundamental fields can be deter- 
mined only by the mass ratios in Eqs. H28|) . For example, 
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without using rn„ and mg, the four mass ratios lead to 



Ti 



rrit I \ vfis I \m 



0.12, (31a) 



D^^['^] (!^) = 0.098, (31b) 



— l^] =0.11 

mb J \mr 



=0.24, 



(31c) 
(31d) 



Here we have used fermion masses at the unification scale 
evaluated in Ref . |35| . From this naive estimation we can 
say that all these VEVs are typically of order 0.1. 

The rest two mass ratios give us the VEV of $ and 
a non-trivial mass relation between up-type quarks and 
charged leptons: 



(32) 



Clearly this relation receives a large correction from un- 
known coefficients in a very complicated way. Addition- 
ally, there are substantial errors in ruu and mt at the 
unification scale. However, we may present the charm 
mass (with the highest power in the relation) avoiding 
such uncertainties, rUc ~ 1.1 GeV, which is consistent 
with the observation within a factor ~ 3. The VEV of $ 
can be estimated by using m„ 



rriu 
mt 



mc 
mt 



or by using me 



= 8.5 X 10" 



1.1 X 10"^. 



(33) 



(34) 



We see that these rough estimates give $ with a large 
uncertainty. 

We have determined all the VEVs from the fermion 
mass ratios, and hence the CKM matrix elements can be 
obtained as predictions. Without the uncertainty of $, 
we can evaluate \Vcb\ 



\Vab\ ^ = 2.7 X 10- 

mb 



(35) 



which reproduces the experimental data [35| with an er- 
ror sm; 
tained 



ror smaller than a factor of two. However, is ob- 



with a large uncertainty due to the range of $ in Eqs. H33|l 
and H34|) . The observed value lies indeed in this range. 
The factorization relation 1)30(1 lead to \Vub\ in the range 



\Vub 



(0.0014-0.011) . 



(37) 



Using the VEV of $ we can evaluate the effective /n- 
parameter through Eq. (|10|l 



rUu 
mt 

mc 
mt 



mc 
mt 

me 
m-r 



( via. 

\mr 



TeV 



0.54 GeV 



(38) 



Although the sixth power of $ in the expression in 
Eq. (|10|l enhances the uncertainly of $, the /x-parameter 
of the weak scale indeed lies in this range. This means 
that the breaking of the symmetry may potentially 
generate the desired value for /i and the charged fermion 
mass hierarchies at the same time. 

To see this point more clearly, we take the eff'ective /i- 
term as an input parameter, e.g., fi = 100 GeV. Then, 
avoiding a large uncertainty $ becomes from Eq. (|31b|l : 



1/6 



= 2.7 X 10" 



(39) 



In this case the mass ratios for up-quark and electron are 
estimated as 



mm 

mt 

me 

mr 



2.6 X 10"*^, 
6.6 X 10"-*. 



which agree with the observed values 



mt 
me 



obs. 



.1 X 10"^ 



= 2.8 X 10" 



obs. 



(40) 
(41) 

(42) 
(43) 



within about a factor of three. Notice that other mass 
ratios in Eqs. (|28|l are independent on We are also 
able to estimate the CKM elements 



\Vus 



0.13, iKfcl- 0.027, iKbl 0.0034. (44) 



It is clear that we obtain rather good predictions of the 
quark mixing angles, which are consistent with the cur- 
rent observation [33 within a factor of two. Therefore, 
the /Lt-parameter of the weak scale and the hierarchies 
of charged fcrmions both suggest the breaking scale 
given in Eq. ((S^ . 



IK. 



m„mt2 mgmb 

1 1 2~ 

mc^ Tot- 4 TOs 

TOt mem^^ mamb _ 



0.40 



0.054 



TO, 



c Tot- 2 



(36) 



B. Neutrinos 

Next, we turn to discuss neutrino masses and mixing. 
We begin with the Majorana masses for the right-handed 
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neutrinos. For example, we take here the VEV of $ to be 
2.7 X 10-3 by taking ^ = 100 GeV, as shown in Eq. 
It is then found from Eq. 1)23(1 that the Majorana masses 
are evaluated as 



M2.3 



$3 = 5.0 X 10^° GeV, 
- $ = 6.7 X 10^^ GeV 



(45) 
(46) 



With these superheavy masses the seesaw mechanism 
works naturally and ensures the smallness of the effec- 
tive Majorana masses for the light neutrinos which are 
almost left-handed states "s^l . The mass matrix for these 
left-handed neutrinos is roughly estimated as 



Hi 





p 


p 


p 


1 


1 


\ p 


1 


1 



where we defined 



(47) 



(48) 



Combining Eqs. I|31b|l . I|31d|) and H39|) the parameter 
p turns out to be 0.85, i.e., at most of order one. 
The obtained result realizes the requirement imposed in 
Eq. 

The particular form of the matrix M^, given in Eq. H47() 
leads to interesting consequences in neutrino properties.^ 
First, we introduce the typical scale for neutrino masses 
nil, which is given by 



Hi 



4.6 X 10^3 eV. 



(49) 



The numerical estimations, which will be performed be- 
low, show that in our model the effective neutrino mass 
matrix M^, generates the neutrino masses (i = 1, 2, 3) 
with a small hierarchy 7713 > 7712 > mi. We may then 
identify the mass squared differences indicated by the at- 
mospheric and solar neutrino oscillations with 5m1^^ = 



and i5mE 



sol 



7712 — mf. Furthermore, the 
typical scale rrii, will be found to approximate the sec- 
ond mass 777.2 ■ The recent analysis [s^l tells us that 
7772 — v^^soi — (7.2-9.9) X 10~3 eV, which is consistent 
with the result in Eq. 1491) within a factor of two. This 
result gives us an important bridge between the neutrino 
masses and the electroweak scale /i-parameter. In fact, 
we can find from Eq. (|10|l with N = 7 that 



P 



(50) 



Thus we can say that if the neutrino masses were smaller 
or lager by an order of magnitude than the observed val- 
ues, it would be impossible to get the desired values for 



The neutrino mass mat rice s in the similar form have been dis- 
cussed in the literature l36l. 



the /i-parameter. Therefore, the neutrino mass scale in 
the solar neutrino oscillation suggests a strong hint for 
the ^-parameter of the electroweak scale in our model. 

The second consequence from Eq. (|47|) is that we have a 
large i'^-Vt mixing angle solution which is required by the 
atmospheric neutrino experiments. This is a direct con- 
sequence of the fact that 52 and 53 have the same charges 
under the gauge group G as well as the discrete symme- 
tries in order to cancel the gauge anomalies and also to 
reproduce the mass hierarchies in down-type quarks and 
charged leptons. Therefore, the model realizes the lop- 
sided structure [27II2II29II30II . 

Finally, we may expect that a large mixing angle for 
I'e-t'/i oscillation as well as a large Ues in the MNS matrix 
since p — 0{1). In fact, one might worry that our model 
would give the same result as in the so-called "anarchy" 
model for neutrino masses [s^ 13^, since the obtained 
mass matrix in Eq. 1)47(1 with p = 1 is identical to that in 
the anarchy model. We find, however, the following dif- 
ferences: (j) The hypothesis in the anarchy model, the re- 
quirement of basis- independence in neutrino flavor space, 
is broken in our model, and it holds only between and 
h'r , which ensures a large mixing in the atmospheric neu- 
trino oscillations, (ii) Even ii p = 1, there are differences 
in the Dirac neutrino mass matrix. Namely, as shown 
in Eq. ((25(1 . we have a structure from both right- and 
left-handed neutrinos. The former one gives a negligible 
effect, but the latter one is important. Due to the Zj 
charges of 5* the 2-1 and 3-1 elements in Mjy are almost 
zero. Only with these missing elements, there appears a 
non-trivial structure in M^. As we will see below, the 
numerical estimates show the clear difference from the 
anarchy model, and in particular there exists the small 
hierarchy in neutrino mixing angles. 



C. Numerical Estimate 

As we have shown above, our mass matrices can pro- 
duce phenomenologically acceptable masses and mixing 
angles of quarks and leptons including neutrinos with 
some uncertainty from unknown coefficients of order one. 
Here we will make a more quantitative estimate by per- 
forming the numerical calculations. 

First, we include the effects of the RGE (renormal- 
ization group equation) evolution. For simplicity, we 
take the boundary at high energy as $ where we define 
all the Yukawa couplings. Then, we solve the one-loop 
RGEs down to the scale Mz and compare with the ob- 
servational data. The threshold corrections from super- 
symmetric particles are included by the universal scale 
MsusY = 1 TeV. 

We treat unknown coefficients of order one as follows: 
Our mass matrices contain 20 independent coefficients 
which are in general complex numbers. Note that some 
are related of each other due to SU(5) groups. These 
coefficients are created randomly in three different ways 
(denoted by the set A, B and C). In all sets the phase 
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TABLE II; VEVs of bifundamental fields and tan /3 obtained 
by minimizing in Eq. I53II . We also show $ and p evaluated 
from these VEVs. 





Ti Di_ T2 D2 tan/? 


$ p 


x' 


Set A 
Set B 
Set C 


0.188 0.103 0.0883 0.284 23.6 
0.175 0.104 0.0836 0.278 23.8 
0.159 0.105 0.0768 0.267 24.0 


0.00271 1.11 
0.00271 1.10 
0.00270 1.09 


0.535 
0.464 
0.404 



of each coefficient is created randomly such that it is 
distributed uniformly in the Hnear scale within the range 
[0, 27r]. The absolute value of each coefficient is created 
also randomly such that it is distributed uniformly in the 
logarithmical scale within the range [1/1.5,1.5] (set A), 
[1/2,2] (set B) and [1/3,3] (set C). 

The fermion mass matrices in our model are deter- 
mined by the six parameters as shown in Eq. (|26|l . In 
this numerical study we set ^ = 100 GeV (at the ultra- 
violet boundary), for simplicity. There are still five free 
parameters - Ti, Di, T2, D2 and tan/? - which have to 
be evaluated from the observational data. Here we use 
the masses for charged fermions m/ and the CKM ma- 
trix elements Vaji, and then the neutrino properties are 
obtained as the outcome. 

We create randomly N — 10^ sets of the unknown 
coefficients, and we calculate m\ and V^p for each set 



and take the mean values in the logarithmical scale as 
follows: 



Vap = exp 



(51) 



(52) 



which should be compared with experimental data. As 
we will present later, some observables such as the CKM 
matrix elements are distributed by orders of magnitude, 
and the mean values in the log-scale are adequate for our 
analysis. 

The five free parameters are fixed by minimizing the 
function defined by 



E 

Oi<OJ 



In 



qj_ 



Of>0 + 



In 



qj_ 



,(53) 



where O/ denotes my or Vap, and and the up- 
per and lower limits of the observable with the 3(t error, 
respectively. In the present analysis we use the masses 
at the scale Mz for charged fermions (except for top- 
quark) estimated in Ref . '35' , and the pole mass for top- 
quark 4Qj. Further, |Ks|, \Vcb\ and |Kb| in Ref. 40] are 
used for the comparison. 

The obtained values for the VEVs and tan (3 are shown 
in Table HTl for the three different sets A, B and C. In Ta- 
ble mil we also show the mean values of fermion masses 
and mixing angles when takes its minimal value. We 
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0.3 
0.2 
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0.5 
0.4 
0.3 
0.2 
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0.3 
0.2 
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FIG. 1: Distributions of charged fermion masses (top) and the 
CKM matrix elements (bottom). Distribution is divided by 
3 only for top-quark mass. The observational data with the 
3(T error are also shown. We use the set B for the unknown 
coefficients. 



see the results do not depend much on the choice of the 
sets for the unknown coefficients. Henceforth, we will 
only use the set B for our discussions. Note that the 
VEVs of bifundamental fields from this numerical esti- 
mate are almost consistent with those shown in Eqs. (|31|l . 
We find tan/3~ 24. 

In Figs, n we show distributions of the charged fermion 
masses and the CKM matrix elements by the N = 10^ 
trials of creating the coefficients. Here only the shape of 
the distribution should be considered, since the hight of 
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TABLE III: Mean values of fermion masses and mixing angles. We also show the observational data with la error (3a" error) 
for charged fermion masses and the CKM elements (neutrino mixing angles). 





Observation 


Set A 


Set B 


bet U 


ruu [MeVJ 


9 00+0.42 
Z.OO_Q 45 


0.859 




0.857 




0.858 


nic [MeVJ 


£?'7'7~l"56 

0''-61 


1.06 




1.05 




1 n^ 
1.0/ 




1 ^7/1 _1_ C C 


181 




179 




1/4 


rud [MeV] 


/I cn+0.60 


3.71 




3.25 




2.69 


rus [MeV] 


no /(+11.8 


51.3 




51.2 




50.2 


rub [CjeVJ 


3.00 ± 0.11 


2.68 




2.75 




2.91 


me [MeVJ 


n /I CT 
0.48/ 


0.832 




0.802 




0. /60 


[MeVJ 


1 no 
103 


72.2 




74.4 




^ A n 


rrir [vjeVJ 


1 '7K 


1.53 




1.57 




l.DO 


n/ 1 

j Vus\ 


n oon ^ n nnoc; 
U.zzU ± 0.002D 


0.187 




0.187 




0.18/ 


IT/ 1 

\Vcb\ 


1 A io^nic\ \y in — 2 
(4.io ± U.iOj X lU 


3.24 X 10 


2 


3.26 X 10" 


2 


Q OQ in — 2 
6.66 X lU 


IT/ 1 

\Vub\ 


(3.0/ ± 0.4/ j X iU 


3.47 X 10" 


3 


3.37 X 10" 


3 


o n'7 V/ 1 n — 3 
3.0/ X lU 


tan^ 023 


0.49-2.2 


1.04 




1.04 




1.05 


tan^ 012 


0.29-0.64 


0.408 




0.350 




0.280 


tan [713 




0.118 




0.0954 




0.0608 


mi [eV] 




4.02 X 10" 


4 


3.63 X 10" 


4 


2.61 X 10"" 


m2 [eV] 




3.11 X 10" 


3 


3.09 X 10" 


3 


3.08 X 10"^ 


ms [eV] 




1.28 X 10" 


2 


1.39 X 10" 


2 


1.69 X 10"^ 


Ml [GeV] 




0.911 X 10 


11 


0.965 X 10 


11 


1.14 X 10" 


M2 [GeV] 




4.36 X 10^'' 


4.84 X 10^^ 


5.34 X 10^^ 


Ms [GeV] 




1.23 X lO^'' 


1.30 X lO^'' 


1.47 X 10^® 



the histogram reflects just the bin size of the horizontal 
axis. The spread of these distributions of masses and 
mixing angles originates in our ignorance of the hidden 
physics to determine the unknown coefficients. 

We can see that our model predicts - of order of mag- 
nitude wise - all masses of charged fermions as well as 
all quark mixing angles by using only five parameters 
(apart from /i = 100 GeV). One may worry that the 
model induces the wrong values for masses and mixing 
angles with some possibility, since they are obtained with 
the distributions, e.g., \Vus \ distributed by two orders of 
magnitude. However, the peaks of the distributions tell 
that the observed values are preferred in the model. In 
this sense, the model can reproduce the charged fermion 
mass spectra and flavor mixing rather well. 

Now let us discuss the neutrino masses and mixing. 
These quantities are obtained as the prediction, since 
we have determined all free parameters from charged 
fermion properties. The distributions of left-handed neu- 
trino masses are shown in Fig. [2| It is found that 
there exists the small hierarchy in neutrino masses (com- 
pared with charged fermions). This feature comes from 
the same reason as in the anarchy model, i.e., the mul- 
tiplication of three mass matrices in the seesaw formula 
(|24|l makes eigenvalues scattered and it is possible to have 
(jnzl'rri'i) ~ 10 [sIjIs^. The peaks of the distributions 
are located at mi ~ 5 x lO"** eV, m2 ~ 3 x 10~^ eV 



and 7713 — 1 X 10^^ eV.^ This result shows that the most 
preferred value of m2 corresponds to the typical neutrino 
mass scale m^, (|49|l : 

= = 4.6 X 10^3 eV - m2 , (54) 

where we have used the fitted values for $ and tan /3 in 
Table HH 

From the small hierarchy in the neutrino masses nn 
the mass squared differences in the atmospheric and solar 
neutrino oscillations are given by 

\j5ml^^ = (3.7-5.8) x 10"^ eV ~ mg , (55) 
^JSn^y = (7.2-9.9) x 10"^ eV ~ ma . (56) 

We see that the most preferred values for and m,2 are 
consistent with the observation within a factor of three or 
so. With this respect, our model is successful to explain 
the observed mass spectra of neutrinos.^ 



We observe numerically that the peak locations of mi and 
(but not m2) depend slightly on the choice of the unknown coef- 
ficients. When we take the wider range of their norm, the peaks 
of mi and 7713 are shifted to the smaller and larger values, re- 
spectively. The diff'erences, however, are within a factor of two. 
Although the factors of three or so are beyond our approximation 
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FIG. 2: Distributions of left-handed neutrino masses rrii (top) 
and neutrino mixing angles tan'^ 9ij (bottom). We use the set 
B for the unknown coefficients. 



In Fig. 121 we show also the distributions of neutrino 
mixing angles, tan^^ij. Although the distributions are 
spread by a few orders of magnitude, we may find char- 
acteristic features in our model. First, similar to the neu- 
trino masses, we find a small hierarchy in the neutrino 
mixing angles, tan^ 023 ^ tan^ 6*12 ^ tan'^6'i3, compared 
with the CKM angels. We should stress that this hierar- 
chy is obtained even when p — 1.1 [see Eq. (j47|l ]. This 
is because of the missing elements in the Dirac neutrino 
mass matrix in Eq. (|25|l . 

The peak of the tan^ 623 distribution suggests that it 
is the maximal mixing which is confirmed in the atmo- 
spheric neutrino oscillations. "'^^ This is the direct con- 
sequence of the fact that the matter fields, and Sg, 
carry the same charges for the product group G as well 
as the discrete symmetry Z7, as mentioned before. The 
peak location of the tan^ ^12 distribution is found to 
be tan^ 6*12 = 0.35, and hence it prefers large values, 
rather than the maximal angle. This is completely dif- 
ferent from the anarchy model where the distributions of 
tan^ 623 and tan^ 6'i2 are the same. 



in this analysis, we may correct them just by changing the overall 
normalization of M^; or Mpf with a factor of \/3 or 1/3, which 
leaves other results unchanged. 

It is found numerically that the peak locations of tan^ 612 and 
tan-^ 613 (but not tan^ 623) depend slightly on the choice of the 
unknown coefficients. The wider range for the norm of the co- 
efficients makes the peaks of tan^ 612 and tan^ 813 smaller. The 
differences, however, are within a factor of two. 



1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 1 1 1 
Set A — 




SetB — 




'"L, SetC - 


=r'i 


'r^fetd^::-!-.. 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

lUeal 

FIG. 3: Distributions of \Ue3\ in the neutrino mixing matrix. 
We use the set A, B and C for the unknown coefficients. 



Moreover, it is interesting to mention that in our model 
tan^ 6*13 < tan^ 9i2 is preferred, and the peak location of 
the tan^ 6'i3 distribution is found to be tan^ ^13 = 0.1, 
which is consistent with the current limit tan^ 6113 < 5.7x 
10~^ jsT] within a factor of two. To be specific, we show 
in Fig.Othc distributions of the \ Ue3 \ element in the MNS 
matrix for the three different sets, A, B and C. The most 
preferred value for \Ue3\ is found to be 0.1-0.3 depending 
on the set. Note that these distributions are plotted in 
the linear horizontal axis, and then we may easily see 
the difference between the three sets. These results show 
that the forthcoming experiments on | Ue3 \ are crucial for 
our model. 

To summarize, our model describes well the masses of 
quarks and charged leptons and also the CKM mixing 
angles when p — 100 GeV. Inversely, the observed mass 
spectra require the /i-parameter around the weak scale. 
In addition, the model predicts the desired features of 
neutrino properties, the correct mass scale and the small 
hierarchies in both masses and mixing angles. We should 
stress again that the typical neutrino mass scale (the so- 
lar neutrino mass scale) also indicates the weak scale fi 
independently on the charged fermion sector. 



V. LEPTOGENESIS 

As studied above, we have determined all the Yukawa 
couplings (or mass matrices) of quarks and leptons in- 
cluding neutrinos, and hence we may calculate various 
quantities in the flavor physics. As an example, we dis- 
cuss here the implication of the model to the leptogenesis 
mechanism |4]| for the cosmic baryon asymmetry. 

The non-equilibrium decays of right-handed Majorana 
neutrinos can generate the lepton asymmetry, if the CP 
symmetry is broken in neutrino sector, because their 
mass terms break the lepton number. The CP asymme- 
try in the decay of right-handed neutrino (here we con- 
sider only the lightest one) is parameterized by ei, which 
is calculated from the interference between the tree and 
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FIG. 4: Distributions of Majorana masses for right-handed 
neutrinos, Afi (top), M2 and M3 (bottom). We use the set B 
for the unknown coefficients. 



one-loop amplitudes for the decay processes [42, |43, |44 
1 

1=2,3 



£1 



Y.lnA{{Y,Y})uY]nx,), (57) 



87r(i;rj)ii 

where Xi = Mf /Ml and /(x;) is defined by 



-^/xih\ ( 1 -I- 



2y^x7 

X,: - 1 



(58) 



Here and hereafter we work in the base of the Majorana 
mass matrix, Eq. H23|) . for right-handed neutrinos being 
diagonal. The distributions of the Majorana masses are 
displayed in Fig.0] We find that typical masses are Mi ~ 
10" GeV, M2 ~ 7 X 10^5 GeV and M3 ~ 1.3 x lO^^ GeV, 
i.e., Ml < M2 < M3. 

Since we have already obtained the Dirac Yukawa cou- 
plings for neutrinos, we are able to calculate the CP 
asymmetry parameter ei. We should note that since all 
the phases of the unknown coefficients are distributed 
randomly in the range of [0, 27r], the sign of ei parameter 
can be positive or negative. In this analysis we neglect 
the sign of ei and discuss only its absolute value. The dis- 
tribution of |ei| at the scale of Mi can be found in Fig.|Sl 
|ei| is distributed by many orders of magnitude, however, 
it is seen that the most preferred value is |ei| ~ 2 x 10^^. 

This lepton asymmetry is produced when the cos- 
mic temperature is T ~ Mi and it is partially con- 
verted into the baryon asymmetry through the elec- 
troweak sphaleron |45| . The baryon-to-entropy ratio of 




mi [eV] 

FIG. 5: Distributions of the CP asymmetry parameter |ei| for 
the lightest right-handed neutrino (top) and the effective neu- 
trino mass mi (bottom). We use the set B for the unknown 
coefficients. 



the present universe is given by 0, 03 

1.5 X 10^3 K|ei| , (59) 



s 



where k is the efficiency factor, which should be esti- 
mated by solving the Boltzmann equations. This factor 
depends on the lightest Majorana mass Mi and the ef- 
fective neutrino mass mi which is defined by [4^ 



TOl 



{MdMI)i 
Ml 



(60) 



Note that M]j is the Dirac neutrino mass matrix in the 
basis where right-handed neutrino mass matrix is diago- 
nal. 

Fig-Elshows also the distribution of mi in our model, 
and we find mi = 0(10^^) eV. According to Ref. 
when Ml = 10^^ GeV and mi = 2 x 10~^ eV, k becomes 
about 10~^. Therefore, the typical values |ei| = 2 x 10^^ 
and K = 10^^ induce the baryon asymmetry 



s 



3 X 10' 



(61) 



which is consistent with the recent observation 
("B/s)|obs. = (9.0 ± 0.4) X 10"" 113 within a factor 
three or so. However, the discrepancy of a factor three 
is of course beyond our estimation of order of magnitude 
wise, and also there is a sizable theoretical uncertainty 
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in the efficiency factor | M l47 | ." 

Tliis is a very encouraging result of the model, i.e., 
the observed cosmic baryon asymmetry is naturally ex- 
plained (with some uncertainty from the unknown coeffi- 
cients) by invoking the leptogenesis mechanism. Unfortu- 
nately, the sign of the asymmetry is not determined in our 
approach. To realize this successful scenario, the maxi- 
mal temperature of the universe should be higher than 
Ml ~ 10^^ GeV so that the lightest right-handed neutri- 
nos were thermalized. Such high temperatures may lead 
to the cosmological gravitino problem in the interes ting 
region of the gravitino mass (of the weak scale) |4a. l49l| . 
To avoid this difficulty we may go into the heavy gravitino 
masses of 0(100) TeV suggested from the anomaly me- 
diated supersymmetry breaking mechanism '50(|.^^ (See 
the discussion in Ref. |12-) 



VI. CONCLUSIONS 

We have investigated in this article a supersymmetric 
model based on the gauge group SU(5)xSU(5)xSU(5) 
with the discrete symmetries Z7 x Z2. The gauge sym- 
metry is broken into the standard model group only by 
using the bifundamental Higgs fields without introduc- 
ing the adjoint or higher- dimensional Higgs fields. This 
point gives us a promising connection to the string theo- 
ries. In fact, all the fields in the model are available in the 
string spectrum with the affine level one. The mass split- 
ting between the weak-doublets and the color-triplets of 
Higgs fields is ensured by the Zj symmetry. The effective 
/i-parameter is then generated by its breaking together 
with the gauge symmetry. 

We have constructed the model for realistic masses 
and mixing of quarks and Icptons including neutrinos. 
The hierarchical structures in the Yukawa couplings orig- 
inate in the breaking pattern of the gauge and discrete 
symmetries. The fermion mass matrices are determined 
by only the six parameters - Ti, Di, T2, D2, $ and 
tan/3 = Hu/Hd - which are all the VEVs of the Higgs 
fields introduced in the model. All the masses and mix- 
ing angles are predicted being consistent with the ob- 
servation including the uncertainty from the unknown 
coefficients of order unity, say, within a factor of two or 
three. 

There have been obtained several interesting features: 
First, the bottom-tau unification, rub — m-r, is real- 
ized within a good accuracy while avoiding the unwanted 
SU(5) mass relations, nig = and md = rrie. Sec- 
ond, right-handed neutrinos carry non-trivial charges of 



The correct value for the baryon asymmetry may be obtained by 
changing the overall normalization of Mo (or Yo) with a factor 
of \/3, which also corrects the discrepancy in the light neutrino 
masses as m entioned in the footnote 1^ 

In Ref. IfSlI a possibility of having high temperatures has been 
discussed. 



the discrete Zt symmetry for the doublet-triplet split- 
ting, and its breaking generates the superheavy Majorana 
masses. Thus, the model offers the natural framework 
for the seesaw mechanism. Third, the Zj breaking gives 
us a non-trivial bridge between three independent facts, 
(i) the /i-parameter, (ii) the hierarchies in the charged 
fermion masses and the quark- mixing angles, and (iii) 
the typical neutrino mass ttLi, (which corresponds to the 
solar neutrino mass scale). It has been shown that both 
of the observational data on (ii) and (iii) point toward 
the Zr breaking scale $ ~ 2.7 x IQ-^ = 6.6 x lO^^ GeV, 
which suggests interestingly the /i-parameter of order of 
the weak scale. Finally, the particular structure in the 
Dirac neutrino mass matrix H25|l generates the small hier- 
archy in the neutrino mixing angles tan^ 023 ^ tan^ 9i2 ^ 
tan^ 013 . The most preferred values of the mixing angles 
indicate that the atmospheric and solar neutrino-mixing 
are maximal and large, respectively. The element \Ue3\ 
in the MNS matrix is typically |J7e3| ~ 0.1-0.3, which 
means that the present model will be tested by the fu- 
ture experiments on Ues- 

We have also discussed the implication of the model to 
the baryon asymmetry of the universe, i.e., the leptogen- 
esis mechanism by the decay of the lightest right-handed 
neutrino which is produced thermally in the early uni- 
verse. Since the model has fixed all the couplings and 
masses of neutrino (including the Majorana masses for 
right-handed neutrinos), we may calculate of order of 
magnitude the baryon asymmetry without any further 
assumptions. Due to the random phases of the unknown 
coefficients the sign of the asymmetry cannot be deter- 
mined. It has been shown that the prediction of the 
model (obtained as the most preferred value) agrees with 
the observation within a factor three or so, and hence 
the model naturally accounts for the baryon asymme- 
try of the present universe. For this successful scenario 
the highest temperature of the universe T > Mi ~ 10^^ 
GeV is required, and the cosmological gravitino problem 
should be avoided somehow. 

Before closing this article, we would like to give some 
comments: It has been found from the fermion masses 
that the gauge symmetry breaking scales, namely the 
VEVs of the bifundamental fields, are 0.08-0.3 ~ (2- 
5) X 10^^ GeV, which is about one order higher than 
the unification scale ~ 2 x 10^^ GeV. This might arise 
a trouble in the gauge coupling unification, although it 
holds at the tree-level. However, we have to take into ac- 
count the threshold corrections from superheavy particles 
beyond the minimal supersymmetric standard model. 
In fact, such particles are indeed present in our model 
and may potentially give rise to the coupling unification. 
Moreover, we keep in mind that there might exist non- 
negligible corrections from the physics beyond our model, 
e.g., the string theories. We should also mention that 
the breaking of our gauge and discrete symmetries gen- 
erates the suppression for the Yukawa couplings for the 
color-triplet Higgs fields and also for the supersymmety- 
breaking masses for the scalar particles. Therefore, our 
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model would avoid the rapid proton decay and also the Acknowledgments 
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